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In the framework of the covariant Bethe-Salpeter formalism we calculate the amplitude of elec- 
trodisintegration of the deuteron near the threshold energy. Due to the small relative momentum 
of the final np-pair, the deuteron break-up process may be sufficiently described by a transition to 
' the dominant So channel into the final state. Therefore only one unknown structure function (form 

^\ , factor) exists. The explicit formula for this form factor is given in terms of partial Bethe-Salpeter 

0^ ' amplitudes. We compare our result to the usual non-relativistic calculations. The essential role of 

P-states as well as their possible link to mesonic exchange currents is discussed. 
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I. INTRODUCTION 



Disintegration of the deuteron (either by photons or electrons) has been and still is a rich source of information 
on the structure of electromagnetic current operators and nuclear dynamics. In the context of photodisintegration a 
recent overview on experimental issues and theoretical methods is given in ref. [Q . 

Here we touch upon the problem of deuteron electrodisintegration near the threshold energy. Due to the small 
■ value of the relative momentum (p* , in the rest system) of the final np-pair it is reasonable to suppose that the 1 Sq 
\ channel dominates in the final state [^J. So the experiment is well suited to study the deuteron wave function in a 
pure way. Recent experiments have been done in Saclay for E np = — 3 MeV 0, and at SLAC with E np = — 10 
MeV Q (E np is the relative energy of the outcoming neuteron proton pair). 

The disintegration process is usually described within a non-relativistic framework including relativistic corrections. 
' Only one investigation is known to have used a relativistic description of the process considered here in the framework 
I | of light front dynamics |5| . 

It is recognized, that the non-relativistic impulse approximation does not reproduce the experimental differential 
cross section even for small momentum transfers. Actually, the minimum that appears in the nonrelativistic impulse 
approximation is not present in the experimental data J3|,|],^]. This minimum, however, is filled with contributions 
of mesonic exchange currents (MEC). In its turn, this filling is a strong confirmation of MEC. Among the MEC 
considered, the most important ones are related to the pion exchange In addition, p-meson exchange currents, 

A-isobar configurations, and final state interactions as well as retardation effects, are to be taking into account dUff). 

Generally, the differential cross section is very sensitive to the choice of the neutron form factor, vertex form factors, 
and the way of constructing mesonic exchange currents. The latter strongly depends on the form of the jNN form 
factor (Fx, or Ge) that influences the results of the cross section ||||. 

The main goal of the present work is to describe the deuteron electrodisintegration in a consistent relativistic 
approach. This will be done with the Bethe-Salpeter equation for the two nucleon system. This way is possible 
to come to general conclusions about the amplitude of the process not seen in the non-relativistic approach. On 
the other hand, the non-relativistic limit will be recovered and some links to the non-relativistic corrections can be 
established. 

In addition, we consider polarization observables in the covariant approach and give some explicit examples for 
polarized electrons on polarized deuteron scattering. The non-relativistic formalism may be found in ref. [ fio| and the 
helicity approach of polarization in ref. fll|| . 

The paper is organized as follows. In the next section the general form of the transition amplitude for deuteron 
break-up near the threshold energy is introduced and differential cross-sections for unpolarizcd and polarized electrons 
and deuterons are given. In section III relations between the asymmetries are deduced. Section IV is devoted to the 
covariant formalism of the Bethe-Salpeter approach, and in section V we evaluate the amplitude in the general form. 
Different approximations are considered in section VI. The connections to nonrelativistic calculations are discussed 
in section VII. The last section contains our conclusions and outlook. 
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II. THE DEUTERON ELECTRODISINTEGRATION NEAR THRESHOLD CROSS SECTION: GENERAL 

FORMULAE 



The general form of the deuteron electrodisintegration amplitude Mfi in the one photon approximation (Fig. [l]) 
may be written in the following way [Q , 



M fi = ie 2 u{k' e ,s' e ) Y l u(k e ,s e ) i (np\j„\DM), 



(1) 



where u(k e , s e ) denotes the free electron spinor with 4-momentum k e and spin s e , and q = k e — k' e is the 4-momentum 
transfer. The hadronic transition matrix element (np\ j^DM) leads from the deuteron state \DAi) with 4-momentum 
K and total angular momentum projection M. to the final np state with 4-momentum P = K + q, where j M is the 
electromagnetic current operator. 

Near the threshold energy the final np state may be assumed to be dominantly in the state |]l3|] and the hadronic 
matrix element may be considered then as a 1 + — ► + transition. Due to the Lorentz transformation properties and 
because of parity and time reversal symmetry, the general form is restricted to only one contribution, i.e. 



(npCSo^j^DM) = ie^ ^ <f ^ V(s,q 2 ) = ^ G m V( Sl q 2 ), 



(2) 



where £^ is the polarization 4-vector of the deuteron, and V(s, q 2 ) a Lorentz scalar real function of invariants 
q 2 and s = P 2 « Am 2 . Below this structure will be shown explicity by straightforward calculations. Note that 
due to the antisymmetric tensor t^ap-y this transition current meets the requirement of gauge independence i.e. 
g At (np( 1 So)\ju\DM) = 0. Thus the evaluation of the amplitude in the framework of the Bethe -Salpeter approach in 
the form of eq. (||) will prove the non trivial property of gauge independence. 

Using eq. (Q) the differential cross-section may be obtained in the standard way see ref. plf . Here it is useful to 
introduce leptonic and hadronic W^ v tensors. Then the differential cross section reads as 



d 2 a 
dE'dVl' 



a 2 I fe e ' I 



(3) 



with a = e 2 /4ir. The leptonic tensor is given by 



2{k e ^k el 



+ 2im e (iuv a m as e- 



The hadronic tensor has the following form: 



= (np( 1 S )\.r\DM)(DM\j^\np( 1 S ))^ 

Using the general form of the hadronic transition current eq. 

= R G» a p aP G* vf} 
where R is a purely kinematical factor. It is given by 

i i IpI 



S(K- 



dk r 



2E p (2tt) 3 2E n (2ir) 3 ' 
the hadronic tensor may be written as 

V(s,q 2 



R = 



8tt 2 2M 



P 



Here m is the nucleon mass and M is the deuteron mass. 
In eq. (^) the density matrix p a p of the deuteron reads 



Paf3 



KgKfj , 

M 2 ' 



1 

2M 



ieapjgK 1 s 5 D 



+-^(-9x1X2 + 



p 
M 2 



)(-9af) + 



M 2 



(4) 



(5) 



(6) 



(7) 



where (W\) a g = itap^xK 1 jM , sd is the spin vector and po is the alignment tensor of the deuteron. Using this 
explicit form of the density matrix the hadronic tensor may then be written as (electron mass is neglected) 
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W$ =\R [g^(<l 2 M 2 - [Kqf) + (K^q v + q M K v )(Kq) - K^K v q 2 - q^M 2 ] V(S, q 2 ) 



Wj$ = \RM ( SD q) ie^ap q° \v(s, q 2 ) 



W^> = R 



K Xl K X2 
M 2 



)[g, u {q 2 M 2 ~{Kq) 2 ) 



+ {K ll q v + q^K v ){Kq) - K^K v q - q^M 2 ] p A D lM V(s, q 2 ) 



The superscripts (u,v,t) denote unpolarized, vector polarized and tensor polarized cases, respectively. 
Thus for the case of unpolarized electrons and deuterons the differential cross section can be written as 

d 2 a (d (7 \ M \p*\m 



dE' e dVL' e 



where 



I^Jm m+K?-™&**i\ Vis,, 2 ) 



( da\ or cos" 

\dnJM = AE 2 sin 4 % 

is the Mott cross section. The normalization conditions of state vectors have been chosen to be for deuteron 

< DP'M'\DPM >= 2E D (2tt) 3 8 mm ,8(P' - P) 

and for nucleons 



(8) 



(9) 



(10) 



< Np, n\Nk, n' >= 2E N (2tt) 3 8^,8{p - k) 



III. ASYMMETRIES FOR POLARIZED DEUTERONS 

With the general form of hadronic tensor W^ u eq. (ph, it is straight forward to calculate asymmetries of the deuteron 
disintegration near threshold, e.g. 

da(tD)-da(UD) 

where da is the differential cross section, f (|) denotes the helicity A e = +1(— 1) of the incoming electron and D the 
polarization state of the deuteron. We assume the initial electron moving along the Z-axis, and 6 e is the electron 
scattering angle. The scattering plane of the electron is in the XZ plane (see Fig^). Then the vectors k e and k' e are 
of the form: 

fc e = (E e ,0,0,E e ), k' e = (E' e ,-E' e sm9 e ,0,E' e cos9 e ). (12) 

First we consider the case of vector polarized deuterons. If the direction of the deuteron polarization is parallel to 
the Z-axis, then 

A _ 3 (E e + E' e )(E e -E' e cos 9 e ) 
11 2 (E e + E' e ) 2 ~ 2E e E' e cos 2 d e /2' [ ' 

where K is the degree of polarization of the deuterons. Note, that the dependence on the form factor V(s,q 2 ) 
disappears. For the case of the backward scattering (Q e — 180°) eq. ( |l3| ) even simplifies to 

A„ = Ik. (14) 

If the polarization of the deuteron is parallel to the X-axis, then 

3 (E e + E>)E> eS m6 e 

± 2 (E e + E' e ) 2 -2E e E' e cos 2 9 e /2' 1 ' 
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These formulae may be generalized to arbitrary polarization direction of the deuteron given by the angles p), 
viz. 

3 (Ee + KWe s in d e sin i3 cos p -\- (E e - E' e cos fl e ) cos g) 
^' ^ - 2 K (W^W^ ' ( j 

Now, consider the case of tensor polarization of the inital target. If the initial deuteron is only aligned due to a 
PDzz component, then the cross section reads 

dcr{pDzz) = da[l + A Z zPD zz }, (17) 
_ 4E 2 + E' e 2 + AE e E' e - AE e E' e cos 6 e + 3E' e 2 cos 29 e 

^-77 



4((E e + E' e ) 2 - 2E e E' e cos 2 e /2) 
where A zz is the tensor analyzing power. For the backward scattering the analyzing power is simplifies 

A zz = 1. (18) 



IV. INITIAL AND FINAL STATES IN THE C O VARIANT FORMALISM 

The main object of the Bcthe Salpctcr approach is the Bcthc Salpctcr amplitude (or wave function), which is 
usually decomposed into a sum of direct products of free Dirac spinors. In order to achieve such a decomposition, an 
additional quantum number for the moving nucleon has been introduced besides spin and isospin, which is the p-spin. 
The operators of spin Si and p-spin Ri (i = 1,2,3) may be represented in the following way: 

Ri =7 5 (k-7 + m)/S fe , 

R 3 =7o(k-7 + m), (19) 



and 



, m ki k • 7 k • 7^ N 

-7075(^7*--^ + — — -=rr), (20) 

y E k E k E k E k (m + E k )' 



where k is the momentum of the nucleon and E k = \Jk 2 + to 2 . The commutation relations for these operators are: 

[Si, Sj] = 2i€ij k S k , [S 2 , Si] = 0, 
[Ri,Rj] — 2ieij k R k , [R 2 ,Ri] = 0, 
[Si, Rj] — 0, 

where S 2 = J2i $i an d R 2 = J2i Ri ■ -"- n the case of k = these operators have a simple form of a direct product of 
matrices acting independently in p and a spaces, viz. 

Si = I p <g> a, R i =pi®I <7 , (21) 

where o~i(pi) are Pauli matrices and Ia(Ip) are unit matrices in the spin (p-spin) space. The eigenvector of operators 
S 2 , S3, i? 2 , i?3 is denoted as E/£(k) with p projection of p-spin and p of the usual spin. Then using 

SsU?(k) = pU?(k), S 2 UP(k)=3UP(k), (22) 
R 3 U^(k) = pE^(k), R 2 U?(k) = 3C/£(k), (23) 

one obtains 

t/;= +1 (k) = up(k), U^ik) = vp(-k), 

where Up(k) and Vp(k) are the free Dirac spinors and spin index p — ±1 corresponds to the spin projection ±5. 
For the Dirac spinors it is convenient to use the following notation: 
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t»^(k) - £(k)u p (0), w„(k) = £(k)« p (0). (24) 
where the boost of particle with mass m and spin-i is given by 

£(k) = / W + " 17 ° , (25) 

with the notation fc = fc M 7 M , and the nucleon 4-momentum k\ = k). In the rest frame of the particles the spinors 
are given by 

«/.(«) = (*n)> «/*(«) = ( ° 



Another useful formula is 

with &2 = (-Efej —ft). The Bethe - Salpeter wave function of the two-nucleon state with angular momentum J and its 
projection M. in the rest frame is written as follows (see also 

*JA<(^,*)=X)fa(*d,|fc|)>%(k), (27) 

a 

where if is the total momentum of the system (in the rest frame K = K^ ' = (M, 0)), and k — (ko, ft) is the relative 
momentum. The decomposition is according to the quantum numbers of relative orbital angular momentum L, total 
spin S, and p-spin collectively denoted by a |l5| . The radial parts of the wave function are denoted by g a (ko, |ft|), 
and the spin-angular momentum part 3% (k) is given by 

3%(k)=t x {Lm L Sm s \JM) (± ml^lSms) 

x(l pi lp 2 \ P m p ) Y LmL (n k ) t/£ (1) (k) U^f\-k), (28) 

where (•]•) denotes the Clebsch-Gordan coefficient jL6l and pi and pi are supposed to be equal ±5. In the following 

formulae, the states with p = 1 are indexed by H — h, or e depending on the quantum numbers of the two particle 

direct product state. The states with zero p-spin are indexed by o respectively. The states denoted by H — or — h do 
not have a definite p-spin. They are connected to the e and o states via 

y + - = -j=(y e + y°), y-+ = -j=(y e -y°). (29) 



The wave function in the form of eqs. (J27-28) is usually referred to as direct product form. However, we find it 
more convenient to work in the matrix scheme and interprete the Dirac index of the first spinor enumerating rows and 
the index of the second spinor enumerating columns. The matrix representation of the Bethe-Salpeter wave functions 
can be obtained from the direct product form eq.(|27j) by transposing the spinor of the second particle. In the rest 
frame of the system this reads for the basis spinors 

^ )(1) (k) t^ )(2) (-k) - c^ (1) (k) E/^f-k), (30) 

which is now a 4 x 4 matrix in the two particle spinor space. So, the two particle Bethe-Salpeter wave function in 
this basis is represented by 

$ JM {K,k) - X JM(K,k) U C = J2 9a(k , |ft|) rfo(k) U c . (31) 

a 

where Uc is the charge conjugation matrix, Uc — 27270 j and is decomposed in the same way as given for 3% 
where eq.([30]) is incorporated into eq.(|2^). Since the matrix notation is not very familiar, we have given an illustration 
of how to calculate the spin-angular momentum part in appendix |a]. 
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To exhibit the p-spin dependence, the spin-angular momentum functions r^(fe) may be replaced in general. Using 
as abbreviation 



fci ± m 

A ± (fc 1 ) = = 32 

y/2E k (m + E k ) 



and the same formula for k\ — > k 2 , one may define 

■ M \ n J = *- M 

where 



r^(fc) = r™(k), (33) 



p a,+4 
1 M 


"(fe) 


= 




pa,— 
1 M 


"(fc) 


= A" 


(h) 


pfi, H — 
1 A4 


"(fe) 


= A^ 


-(h) 


pa, — 1 
1 .M 


"(fc) 


= A" 


(h) 



I±™ f& (fc) A-(fea), 
i±^f^(fc) A+(fcO, 

A(1 2 70) f^(fe) A-(fa), (34) 

with a 6 {L, 5, J}, A = —1 for the 1 Sq channel and A = +1 for the 3 Si — 3 D\ channel. The expressions for the 
functions T% l (k) are given in tables || and 0. 
The adjoint functions are defined through 

r5t<(ft) = 7o [TX,(fe)] t 70, (35) 
with the orthogonality condition given via: 

/ dk Tr{T a M \k)T a MI (y)} = &MM'5aa>- (36) 

In addition, for identical particles the Pauli principle holds, viz. 

XJM (K,k) = -P 12XJ M(K,k) = (-l) I+1 U c \xJM(K,-k)f U c - (37) 

where I denotes the channel isospin. This induces a definite transformation property of the radial functions g a (ko, \k\) 
on replacing fco — > — ko, which is even or odd, independend of a. Moreover, since the 3 0-)ppip3 amplitudes do not 
have a definite symmetry on particle exchange, it is more convenient to introduce rj^, or instead of using 
or rj^f + . These functions do have a definite even(e) or odd(o) p-parity which reflects a definite symmetry behaviour 
under particle exchange. 

Subsequently, we give the explicit forms of the partial amplitude for the 1 S'o and the 3 Si channels. 



A. The 1 S channel 

For the two nucleon system in the 1 5'o channel the relativistic wave function consists of four states, i.e. 1 £ , f+ , 
1 S , 3 Pq, 3 Po, ( 3 Pq , 3 Po + )> enumerated by 1, ... ,4, respectively, in the following. The radial functions of the 
partial wave decomposition of the Bethe-Salpeter amplitude are denoted by g a (k , |fc|) — > 4>i(k , |k|), where 



Introducing Lorentz invariant functions bi(Kk, fc 2 ), the Bethe-Salpeter amplitude xoo(^ 5 fe) is given by (arguments 
of bi suppressed) 

47T xoo(K, k) = foi7 5 + b 2 — (pi75 + 75P 2 ) 

m 

,pi-m p 2 +m pi -m p 2 + m 

+M 75 - 75 ) + 04 75 , (39) 

m m mm 

where p\ y2 = K/2 ± k. The functions bi connected to 4n are given in Appendix [B]. Note, that only b 2 and </>4 are odd 
with respect to fco — > —fco, an d all the other functions are even. 
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B. 3 Si - 3 Di channel 



In the deuteron channel the relativistic wave function consists of eight states, i.e. 3 Si + , 3 S 1 3 Df + 7 3 D 1 , 
3 Pf, 3 Pf, l Pf, ^f, ( 3 F 1 +_ , 3 Pf + , 1 Pf + ), enumerated by 1, ... ,8, respectively, in the following. The radial 

functions of the partial wave decomposition of the Bethe-Salpeter amplitude are denoted by g a (ko, \k\) — > ipi(ko, |k|). 

The covariant form of the Bcthe Salpetcr amplitude xim(K, k) with eight Lorentz invariant functions hi(Kk,k 2 ) 
is given by (arguments of hi suppressed) 

4. X1M (K, k) = h,i M + /> 2 ^ + , 3 ( hzH^ + (40) 

TO y TO TO y 

+/i 4 h /i 5 - 



TO / TO \ TO TO 

. , (§i-P2- 2to\ pi - m / - (^m)\ k + 

+ra 6 1 n-7c;x + «8- 

y to y to to y 

The relation of functions /ij to ipi are given in appendix |b[ 

Note now, that /13, /14 and -05, are odd, and all the other functions are even under k,Q — > —kg. 



V. THE ELECTROMAGNETIC CURRENT OF THE HADRONIC SYSTEM 

The electromagnetic current matrix element (np^ So)\j ^{DM) will be evaluated in the relativistic impulse approx- 
imation, see Fig. ||. Within the framework of the Bethe-Salpeter (BS) approach using the Mandelstam procedure of 
constructing of the electromagnetic current operator , it can be written as 

(np( 1 S )\j li \DM) = i I d 4 kTr(xs(P,P)^Hu)XM(K,k)(^-k + m)\, (41) 



where Xs(P,p) is the BS amplitude of the 1 Sq state of the np-system (similar to the nonrelativistic wave function 
of the continuum), and XM(K,k) is the BS amplitude of the deuteron and p = k + q/2. The vertex of the jNN 
interaction is 

which is chosen to be on the mass shell. The isovector form factors of the nucleon appears due to summation of the 
two nucleons. Instead of Dirac (F^ V \q 2 )) and Pauli (F% (q 2 )) form factors one may alternatively choose Sachs form 
factors. The relation between them is given by 

G B (g 2 )=iM<z 2 ) + £^ 2 (<z 2 ), 

G M (q 2 )^F 1 (q 2 ) + F 2 (q 2 ). (42) 



In order to extract V(s,q 2 ) from the general expression for the 1 + — > + transition from eq. (p4[), we first perform 
the trace [Q. With respect to the Lorentz structure of the eq.(f4l|) seven integrals have to be evaluated, which are 
expressed through Lorentz scalar quantities abreviated by (. . .) (the index M. on the polarization vector £ is suppressed 
for simplicity). 

ix = fVjn J d^k (...), 

I 2 = iW 7 r^ / d A k (...) k\ 
2s = £V / d 4 k (...) k\ 



I 4 = ie^fr J d A k (...) k a k s , 
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X 5 = ie a ^ S Cq P K^ J d A k (...) k 5 k^ 

X 6 = KaPyS Cq P K^K^ J d 4 k (...) k 5 , 

It = *e aPlS Cq P K^q^ J d A k (...) k s . (43) 

From inspection of eq. (|4^) it becomes clear that besides Xi, which needs no further consideration, two generic types 
of integrals have to be evaluated, viz. 

d 4 k (...) k a =C iqa +C 2 K a , (44) 
d 4 k (...) k a kp = D 1 M 2 g af3 + D 2 q a qp + Dz{q a Kp + K a qp) + DiK a Kp, 

where we have already indicated that after integration the expressions should depend on the external 4-momenta only, 
i.e. the transition momentum q and the total deuteron momentum K. Due to the antisymmetric tensor, the number 
of terms is reduced substantially, and all the terms with Z)j, i > 1 vanish. So, besides X\, it is necessary to evaluate 
the following integrals: 



J 2 = ^a/3 7 r^V J d 4 k (...) Cl (k,q,K), 

X 3 = »W7 rY^ 7 / d 4 k (...) c 2 {k, q, K), 



Ii = ^W7 £V# 7 M 2 / d 4 k (...) d(k,q,K), 



X 5 = ie. 



tfSlll Cq 13 ^ M 2 J d 4 k (...) d(k, q, K), (45) 

and Xq = = 0; Ci and d are scalar functions given in Appendix 

Integration and comparison with the structure of the transition matrix element given in eq.(|^) finally leads to 
determine V(s,q 2 ). The result is 

V(S, q 2 ) = Vt (s, q 2 )F[ V) (q 2 ) + V 2 (s, q 2 )F^ V \q 2 ). (46) 
The expressions for V\^ 2 {s, q 2 ) are rather lengthy and therefore shifted to appendix ^|. 



VI. RELATIVISTIC PLANE WAVE IMPULSE APPROXIMATION 

It is instructive to approximate the final np state by a plane wave in order to get a better impression of the terms 
that contribute to the form factor V(s, q 2 ) defined in eq. (^) and evaluated in the last section. We use the covariant 
form of the BS amplitude of the np-system and the deuteron. In the plane wave approximation the l S$ channel radial 
amplitudes in the rest frame of the pair are given by 

0i 5 + (PO, \p\) = —j= ~T2 S (Po) S (\P\ - \P*\)> 
V47T V 



■s-(Po, \p\) = 0, 
»p « (Po, \p\) = 0, 

.j~(po>|)=0, (47) 



where po, \p\ stand for the relative momentum of the pair in the rest frame and \p* \ — \J s/4 — to 2 . 
The calculations have been done in the laboratory system. Writing the (5-functions in eq. (|47|) as 

Pp = _ kq -l,2Mq^l,2q 2 

M + q V M + q h 

y/s/4 - m 2 (M + q ) s - 2E k (M + q ) . 



Si^JI^P- ffiW) = EkMq] S(cosO k agf]^), (48) 
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one can perform integration in eq.(ftq) and set 



M 

ko = k = — - E k , 



s-2E k (M + q ) 
2\k\\q\ 



(49) 



The condition | cos 6u\ < 1 restricts the integration limits of \k\ and the following expression for V\^(s,q 2 ) can be 
obtained 



where 



and 



1 



4tt 2\q\\p* 



\k\ r 



\k\, 



d\k\ L±I h2 ( q 2 , s ,\k\ 



/a a.r .111:1:11 



1 , Q S 2 _u 9 

1 + TVl" m ± T 



y/((M 2 + S - q 2 ) 2 - 4sM 2 ) 

2M 



h( S ,q 2 ,\k\) 



(m + E k ) 
V2a 



[s(s + AEl) - 4mq 2 (m + 2E k ) - 4sE k b] 



9s(K,\k\) 
(M - 2E k ) 



h(s,q 2 ,\k\) = 



2k 2 

•\/3a 
~2M\k\ 
a 



[s(m + 2E k )(s + 4E 2 - 4E k b) + 4mq 2 (E k (E k + 2m) + k 



2y] 9D(ko,\ k \) 
(M - 2E k ) 



2m(m + E k ) 



[s{s 4E 2 ) 4m 2 q 2 ] |fe|) - g° t (k* , |fe|)), 

[s[s(M - 4{E k + m)) - 4E k M(E k + 2m)} + 4m 2 Mq 2 



-is (Mm + 2E k (E k + m))b] 



9s(K, |fcj) 
(M - 2E k ) 



V2a 
4mk' 



■ [s 2 (M(E k + 2m) - 4k 2 ) + 4MsE k {2m(E k + m) - El) 



As(Mm(m + 2E k ) - 2E k k 2 )b + 4Mm 2 q 2 {E k + 2m)] ff 



9d(K, \k\) 



(M - 2E k ) 



V3a 
' M\k : 
V3V2a 
4m\k\ 



(s - 2bE k )(s - Mb)(g e t (k* , \k\) - g<Z(k* 0) \k\)) 

[s(s + 4E 2 - 4bE k ) + 4m 2 q 2 ] |fc|) - g°(k* , \k\ 



(50) 



(51) 



(52) 



(53) 



with a — l/4Msq 2 , b = y/s + q 2 , &q = M/2 — E k , and g a (ko, \k\) are the vertex functions of the partial states of the 
deuteron. The indices S, D denote the 3 S^ , 3 Df states, resp., and index s(t) denotes the singlet (triplet) states and 
o(e) the odd (even) state with respect to the p-spin. 

As it can be seen from eqs. (|52-|53) only the (++) or ( — h) components of the deuteron function are present. This 
can easily be understood by writing the matrix element eq. (^IJ) in the laboratory system and using wave functions 
of the form eq. ([m]). The expression for the inverse propogator in eq. (Ell) is of the following form: 



y - k + m = [(ki - m){— - fc + E k ) + (k 2 + ™)(— - k - E k )] , 

with ki = (E k , k) and k 2 = (E k — k). Thus, for the ( ) component we get 

K 

XM{K,k){— - k + m) ~ 

(k - m)^Y^T M (k)(ki + m) [{ki - m){— -k + E k ) + (k 2 + m){— - k - E k )] = 



(54) 



(55) 
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because of kg — M/2 — E). and (hi + m)(fci — m) — 0. The same is valid for the (H — ) components. 
Then the transition form factor V(s, q 2 ) has a specific simple form in the case of s = Am 2 , viz. 



V(Am 2 ,q 2 ) = 



1 



4?r 1 2Mm 

V3 



1 {F[ V \ q 2 ) + Fi V \ q 2 )) gg(fc "' l9l/2) '^ (fc °' lgl/2) 



M -2 y / m 2 + q 2 /4 



2M 2 \q\ 



(*i (v V 



M 
2m 



Am 2 



FP( q 2 ))(gt(k*,\ q \/2)-g?(K,\ q \/2)) 



(56) 



where fcg = M/2 - y/m 2 + q 2 /A. 



VII. NON-RELATIVISTIC INTERPRETATION AND THE PAIR CURRENT 



Disintegration of the deuteron at threshold has been a cornerstone to establish subnucleonic degrees of freedom 
in the electromagnetic current operator. Compared to the data the non-relativistic impulse approximation fails 
completely, and additional corrections, viz. two body currents have been introduced. On the other hand, relativistic 
impulse approximation partially includes two body currents. Recently, this has been shown explicitly in the context 
of the light front approach to the two nucleon system Q. Here we consider the covariant approach and like to 
offer a motivation why the two body currents are related to the P-wave components which are not present in the 
nonrelativistic treatment. 

As a demonstration we compare the result of the relativistic calculation of the deuteron disintegration amplitude 
for the case p* = with the nonrelativistic one. In this case the nonrelativistic formula is usually parameterized in 
term of the following reduced matrix element 



f So \\T Ma9 \\D} = 



1 



r\q\ 



+ Av)^ 91nn 2 



Ait m\q\ 



d\r\ 



(1 + Mkl) (u s (|r|) + ^(M)) ii(^) 



(57) 



where us and up are the nonrelativistic S— and D— wave functions. The first term in eq. ( p7|) corresponds to the 
nonrelativistic impulse approximation and the second term - to mesonic exchange currents. Within the nonrelativistic 
calculations the nucleon form factor t^ v \q 2 ) in eq. (|57]) has been subjected to discussion. Two possibilities have been 

considered, i.e. ^ v Hq 2 ) — F^ V \q 2 ) || and ^ v \q 2 ) — G^\q 2 ) ||. Within the covariant approach used here, form 
factor ^ v \q 2 ) cannot be chosen freely, that will be explained below. 

In the limit fcj — > it is possible to relate the nonrelativistic wave functions to the relativistic vertex function, 
viz. Si 



0, |fc|) -» a m (2E k - M) u S ,fl(|fe|), 



(58) 



9s,d{K 

with a = Ak\[2M. I m. Inserting eq. 



|) into eq. ( |5q ) one recovers the structure of the nonrelativistic result as given in 
cq. (57). The term proportional to (F\ + F2) is determined by the positive energy components of the Bethe-Salpeter 



amplitudes, g$ and go- We argue that they may be interpreted as the nonrelativistic impulse approximation. Then 
the mesonic exchange current contribution eq. (|57j) may be related to the remaining terms that depend on the P-wave 
components only. 

Since within the covariant approach the sets of form factors {Pi, P2} and {Ge, Gm} are completely equivalent, 
the discussion whether to use Pi or Ge in the form factors of the mesonic exchange part, is obsolete. Our result is 



00 



M 
2m 



Am 2 



M 
2m 



Am 1 



Am 2 



(59) 



This is not the case for nonrelativistic approximations where covariance is broken and higher order corrections in ql m 
are not treated properly. In this case some physics aspects may be described better by using either one of them [§!). 
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With the definition of eq. ( (42| ) and using Ge we may express the form factor dependent part of the second term 
in eq. ( |56|) belonging to the mesonic exchange currents. For simplicity we neglect the deuteron binding energy (viz. 
M = 2m), and 0(q 4 ) terms in the expansion of the square root and write 

W)=ifV) + ^^V)+o(<? 4 ) 

= G ( P-£-A V) (l 2 ) + 0( q 4 ). (60) 

VIII. CONCLUSION AND OUTLOOK 

Although the nonrelativistic limit of the Bethe-Salpeter equation does not exist in the general case, it is possible 
to connect relativistic and nonrelativistic calculations of the deuteron disintegration amplitude. The nonrelativistic 
impulse approximation corresponds to the states with positive energies. Corrections due to the P-components are 
presumably connected to mesonic exchange currents. In further communications we shall give an analytic proof that 
the main contribution of the P-states to the amplitude coincides with the non relativistic pair current. This conclusion 
is supported by the calculation of the magnetic moment, where it has been shown that the size and sign of the P-state 
contribution are the same as in the pair current contribution within the nonrelativistic framework. 

It has been shown that the Bethe-Salpeter approach allows to calculate the amplitude of the disintegration process 
in a covariant form. This is a great advantage, since it gives an opportunity to consider some old problems from a new 
point of view. One of the problems is gauge independence of the current matrix element. It has not been solved for 
the non relativistic treatment. Therefore the problem of choosing the correct form factor in meson exchange currents 
stays persistant. A fully relativistic approach is free from these biases. Another interesting result unseen in the non 
relativistic approach is that simple relations between polarization observables can be achieved. 

Usually the Bethe-Sapleter equation is solved after Wick rotation. The formula of the present calculation has 
shown that the dependence of the vertex function on the variable k is sufficient. This is similar to calculations in the 
context of the pd elastic scattering [^3| where it has been motivated that the dependence on ko does not substantially 
influence the results by choosing ko = 0. However, in the case of deuteron disintegration in the relevant regime of 
momentum transfer, q 2 < 3GeV 2 , the values of ko vary from to 1 GeV. The dependence of the vertex function on 
fco is displayed in Figure ||. The difference of the vertex functions with fco = and fco = k* ^ can be as large as 
20%. 
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APPENDIX A: CALULATION OF SPIN-ORBIT DEPENDENCE 



Here is an example how to calculate the spin - angular momentum part of the wave function for the 3 S^ + - state 
where the spectroscopic notation 2S+1 Lj of ref. [|l5| is used. In this case 



47rr£ (k)U c = J2 (|/*ii/*a|lM)u w (fc)u^(-fc) 

= £(fc) (b^2\lM) ( ^ ) ( xl ) C T (-k) 



N „ ( E (|mi|m 2 |im) Xmi x^ o\ 

= £(fc) /ii/12 £ f— «) 

V o o/ 

= t{k) I±2> * ( o - CT n ^ ) ( ^ C T (-k) 

= £(fe)i±^-L(- 7 ^)£(fe)C/ c 

(m + fci) 1 ^ ) 7 ° £m ("i - fca) t/c, 



2E k (m + E k ) ^2 v ' 2 

Here we have used the identity \[2 (2MiiM2|lA4)x Ml X„ 2 = (er • £ M ) («cr 2 ), where fc 2 = (Ek,—k), £ M is the 
polarization vector of the spin-1 composite system with the components in the rest frame given by 

C + i = (-l,-i,0)/V2, €_i = (l,-i,0)A/5, Co = (0,0,1). (Al) 
and the 4- vector Cm — (0, Cm)- 



APPENDIX B: RELATION OF COVARIANT FORM TO PARTIAL WAVE EXPANSION 

The relation between functions hi of the covariant expansion to the functions <pi of the partial wave expansion is 
given by, for the 1 Sq channel: 

V^2 1 1 

bi = — — ai[Df 0i - D^4> 2 ] - -a 1 a 2 D 2 (t)3 + — a 2 k$(j)± (Bl) 
lb o ivi 

b 2 = —t^fa (B2) 

\/2 1 1 

6 3 = — aim 2 [0i - 02 ] - 2 aifl2 ^^ ~ 2m 2 )0 3 + ^ra 2 fc o 04 (B3) 

64 = -— aiTO 2 [0! - 2 ] - aia 2 m 2 4> z (B4) 



where 



a x = 1/ATEfc, 

a 2 = m/|fe|, 
L>± = (M ± 2£ /c ) 2 - 4(4m 2 + /c 2 ), 
L> 2 = M 2 + 12E% - 16m 2 - 4A& 



for the 3 Si — 3 D\ channel: 



hi = —ai[D^^i + D^ip 2 ] - —ai[D^ip 3 + Dfipi] + —-:a 2 k^ 5 - — aia 2 D 2 i/j 6 (B5) 
lb lo IM lb 



hz = ^a 1 a 3 [D 3 0! + D3V2] + y^ a i a 2 [-°4 V>3 + D 4 ip 4 ] - j^a 2 k il; 5 (DC) 
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where 



16 E k 
h3 = --g-a 2 ip5 

/i 4 = — aia 3 mko[ipi - ip 2 ] - -aiaj{E k + 2m)k [tp 3 - -04] ^-aia 2 mkoip 7 + ~~^ a ^ 



V2 



1 



V6 



V6 



/i 5 = — —aim [ipi + ip 2 ] + 2°i TO bPa + V^] + jTTa 2 k ip5 yaia 2 (£ fc - 2m )~ip 6 

he = —aia^mlD^^i + DqiP 2 ] + -aia 2 [D^ tp 3 + Df ^4] — aia 2 m 2 V>6 + —p 7^02^7 

o o Z 4 iife 

/i 7 = ^aim 2 [0i + ^2] - ^ainJ 2 ^ + 04] - ^-aia 2 m 2 ip 6 



/ig = --yaia3«T- 2 [0i + ^2] - ^a 1 alm(2E k + m)[ip 3 + ^4] + ■^-a 1 a 2 m 2 tl> 6 



a 3 = m/(E k + m), 
Df = (M ± 2E k ± 4m) 2 - 4fc 2 , 

= [(M ± 2E k f - Akl\{2E k + m) - 16(m ± M)fc 2 , 
D 5 = Af 2 - 4£ 2 + 16m 2 - 4fc 2 , 
Df = 4m + 2E k ± M, 

= ±M(E k + 2m) + 2 (2m 2 + 2mE fc — i? 2 ). 



APPENDIX C: THE TRANSITION FORM FACTOR 



(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 
(B13) 



The eq. ( f46| ) for V(s, q 2 ) given in section III is 

V(s,<?) = VW)*fV) + V 2 ( S ,q 2 )F^ V \q 2 ). 

The expressions for Vi t2 (s,q 2 ) are 

Vi(8, q 2 ) = | d 4 fc q 2 , k) + V 2 {1) (s, q 2 , k) + ^(s, q 2 , k) + V^( S , q 2 , k)) 

V 2 (s, q 2 ) = j d'k (V/ 2) ( S , q 2 , k) + V™(s, q 2 , k) + V 3 (2) ( S , q 2 ,k) + v} 2 \s, q 2 , k)) 



where 



(i) 



V, 



(1) _ 



V, 



(1) _ 



v. 



(1) 



Aujicibi h 3 + hs + 2/17 

-2 (2(a - c 2 ) + l) h + 8c 2 h 3 - Aui 2 dh 4 + 8 hci - c 2 + lj h 5 - Aui 2 dh 6 
+ y (^2(4fc 2 + 12m 2 + M 2 - 4(Kk))ci + 2(4fc 2 - 4m 2 - A(Kk) + M 2 )c 2 
+Ak 2 + 12m 2 + M 2 - 4(JTjfe)W - 8u 2 dh 8 

tuib 3 -2^2(ci - c 2 ) + ljhi - 8[2ci - c 2 "jh 3 - Auj 2 dh A + 8^1 - c 2 "jh 5 - Au 2 dh 6 
-y (2(M 2 + 12m 2 + 4fc 2 - 4(Kk))a + 2(4m 2 - 4fc 2 - Af 2 + A{Kk))c 2 
—M 2 - 12m 2 - Ak 2 + A{Kk)^jh 7 - 8uj 2 dh 8 

ujibi -4(2( Cl - c 2 ) + l) hi - Lj 3 [(l2m 2 + M 2 - 8(kq) + Ak 2 - A{Kk) + A{Kq))ci 



(CI) 
(C2) 
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(2) 



+2{M 2 - 4m 2 - 4k 2 )c 2 - 4(Kk) + 8fc 2 - 16dM 2 )) h 3 - 8uj 2 dh 4 

+uj 3 ((12m 2 - 4{Kq) + M 2 + 4k 2 + 8(kq) - 4(Kk))a + 4(2(Kk) - 4m 2 - M 2 )c 2 

+M 2 + 12m 2 - 4k 2 + 16dM 2 )/i 5 - 8us 2 dh 6 + 4^2(feg) - (Kqfjci 

+2lo 3 {{4{Kk) - 4m 2 + 4k 2 - 3M 2 )c 2 - 4k 2 + M 2 + 4m 2 + 16dM 2 ) h 7 

-UJ2 (M 2 + 12m 2 + 4k 2 - 4(Kk)^j dh 8 



vl 2) = 



K (2) = 



4 

Wl&2 



(l - 2c 2 )fti + 4c 2 /i 3 - 2w 2 rf/i 4 + 4^1 - c 2 )/i 5 - 2w 2 d/i 6 
(2(4fc 2 - 4m 2 - 4{Kk) + M 2 )c 2 + 4k 2 + 12m 2 + M 2 - 4(Kkfj h 7 - 4uj 2 dh 8 
4c 2 hi - 2uj 2 dh 2 + 2uj 3 ({4(kq) + q 2 )c t + 4c 2 fc 2 - 2fc 2 + 6dM 2 ) h 3 - 4uj 2 dh A 
+2uj 3 [{4{kq) + q 2 ) Cl + {4{Kk) - M 2 - 4m 2 )c 2 - 2fc 2 + 6dM 2 )/i 5 + 4iu 2 dh 6 
+LU 3 (l6 Cl (kq) + 4c iq 2 + (12(Kk) - 4m 2 + 4k 2 - 3M 2 )c 2 - 8k 2 + 24dM 2 ^h 7 
Jm 2 - 4m 2 + 4fc 2 - 4(Kkfj dh s 

4(l- c 2 ) hi - 2u 2 dh 2 + 2u 3 ((2(Kq) + q 2 ) Cl + {M 2 - 4m 2 )c 2 
2{Kk) + 2k 2 - 2dM 2 ) h 3 + 4u 2 dh A + lo 3 (2(2(Kq) + q 2 ) Cl + 2(4m 2 + M 2 )c 2 
12m 2 - M 2 - 4dM 2 ^) h 5 + 12uj 2 dh 6 + (4(q 2 + 2{Kq)) 



LU 2 

2 



)ci 

+ (3M 2 + 4m 2 - 4k 2 + 4{Kk))c 2 - 8m 2 - 2M 2 - 8dM 2 ^jh 7 



+ 



y (28m 2 + M 2 + 4k 2 - 4(Jffc)) dh s 
LUibi [y (s,(2(kq) - (Kq))a + 2(4k 2 - 4m 2 + 4{Kk) - 3M 2 )c 2 
-12k 2 + M 2 + 12m 2 + 4{Kk) + 32dM 2 )/ii - 4u 2 dh 2 + (4{q 2 + (Kq) + 2(kq))a 
+ (M 2 + 4k 2 - 4m 2 + 4{Kk))c 2 - 4{Kk) + 8dM 2 ^jh 3 

-y (M 2 - 4m 2 + 4fc 2 - 4(Kkfjdh 4 - oj 3 (^4(2(kq) - 3{Kq) + 4{kq) - q 2 )a 
+ (4k 2 + 4{Kk) - 4m 2 - 7M 2 )c 2 - 8fc 2 + 2Af 2 + 8m 2 + 4{Kk) + 24dM 2 ) h 5 
+ y (28m 2 + M 2 + 4k 2 - 4(Kk)^jdh 6 + y ({M 2 {Kq) + 12g 2 m 2 + 2M 2 (kq) 

-8{kq){Kk) + q 2 M 2 - 8(kq)m 2 + 28{Kq)m 2 + 4q 2 k 2 + 4{Kq)k 2 + 8k 2 (kq) 
-4(Kq){Kk) - 4q 2 {Kk)) Cl + 2(16m 4 + 16fc 4 + M 4 - 16(i^fc) 2 - 32fc 2 m 2 
+8Af 2 fc 2 + 56M 2 m 2 )c 2 + 96fc 2 m 2 - 8M 2 fc 2 - 64(if fc)m 2 - 40M 2 m 2 - M 4 

+ 16(iffc) 2 - 16fc 4 - 80m 4 + 16(4fc 2 + M 2 - 4(Kk) - 20m 2 )M 2 d) h 7 
+w 5 (m 2 - 4{Kk) + 4k 2 + I2m 2 ^dh s 



uj^ = 1/m 4 , U5 = M 2 /m 4 . 



with the abreviations 

oji = 1/m, uo 2 = M 2 /m 2 , 0^3 = 1/m 2 , 
Functions c, = Cj(fc, q, K) and d = d(k, q, K) are given by 

M 2 (kq) - (Kk)(Kq) 



d(k,q,K) 



M 2 q 2 - (Kq) 2 
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Co(ka K) _ (Kk)q* - (Kq){kq) 

c 2 {k,q,K)- MW _ {Rq y (03) 
_ (iffc)V + M 2 (fcg) 2 + {Kqfe - M 2 g 2 fc 2 - 2(Kq)(Kk)(kq) 
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TABLE I. Spin angular momentum parts Tq for the 1 So channel 
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Vs^ f o 






-75 




3 Po 


-l/2(fci - fc 2 )|fc 
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TABLE II. Spin angular momentum parts for 3 Si — 3 Di channel 
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V8tt T m 




3 s 1 






-71 




6m + §(&i - k 2 )(M M 


)fc- 2 _ 




3 ^ y/i 




|fe|" 


i 


'Pi 







1G 



D: K,M np: P = K + q 

FIG. 1. The one photon approximation. 
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D: K,M np: P = K + q 

FIG. 3. The impulse approximation. 




np: P = K + q 



D: K,M 

FIG. 4. The plane wave approximation. 
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p[GeV] 



FIG. 5. The fco dependence of the vertex functions. The short dashed line is the relativistic vertex function with fco = 0, 
the long dashed line is the relativistic vertex function with fco = M/2 — Ek, the solid line is the nonrelativistic vertex function. 
The functions are fron the Graz II potential [24] 
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